model1 = @(mu,x) 1-mu*x^2;

N_pass=300;
N=800;
x=-0.1;
mu_N = 6e2;
% mu_beg = 0;
% mu_end = 2;
mu_beg = 1.74;
mu_end = 1.86;

res = zeros(N-N_pass,mu_N);
for muIdx = (0:mu_N-1)
    mu = muIdx/(mu_N-1)*(mu_end-mu_beg)+mu_beg;
    for nn=1:N
        x = model1(mu,x);
        if nn>N_pass
            res(nn-N_pass,muIdx+1) = x;
        end
    end
end

fig = figure('Visible','off');
hold on;
for muIdx = (0:mu_N-1)
    mu = muIdx/(mu_N-1)*(mu_end-mu_beg)+mu_beg;
    mu = mu*ones(N-N_pass,1);
    scatter(mu,res(:,muIdx+1),1);
end
title('chongkou model1');
xlabel('\mu');
ylabel('x');
grid on;
exportgraphics(fig,'test1.jpg');
close(fig);